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ABSTRACT 



The evolution of the magnetic field of an accreting magnetic white dwarf with 
an initially dipolar field at the surface has been studied for non-spherical accretion 
under simplifying assumptions. Accretion on to the polar regions tends to advect the 
field toward the stellar equator which is then buried. This tendency is countered by 
Ohmic diffusion and magneto- hydrodynamic instabilities. It is argued that if matter is 
O ' accreted at a rate of Merit ~ 10^^ g s~^ and the total mass accreted exceeds a critical 

value AMcrit ~ 0.1 — 0.2Mo, the field may be expected to be restructured, and 
^ ■ the polar field to be reduced reaching a minimum value of ~ 10'^ G (the "bottom 

' field") independently of the initial field strength. Below this critical accretion rate, 

the field diffuses faster than it can be advected, and accretion has little effect on field 
strength and structure. 

' In polars, where the magnetic field strength 10^ — 10* G) is strong enough to 

lock the magnetic white dwarf into synchronous rotation with the orbit and a disc does 

(N . 

not form, magnetic braking is severely curtailed as the stellar wind from the secondary 
becomes trapped in the combined magnetosphere of the two stars. The mass transfer 
rate in such systems is typically ~ 10^^ g s""'^, and field restructuring is not expected. 
|0 I In systems with fields not strong enough to achieve synchronism and where accretion 

. occurs via a truncated disc (the intermediate polars) , normal magnetic braking may be 

I expected. The mass transfer rates are then typically ~ 10^^ g s~^ above the 2 — 3 hour 

Cataclysmic Variable period gap, and thus a significant reduction of the polar 
field strength could occur if such a system accumulates the required critical mass 
Merit- However, due to mass loss during nova eruptions, only a small fraction of such 
}2| ■ systems (those that first come into contact at long orbital periods) may accumulate 

' sufficient mass to reach the bottom field configuration. We argue that the observed 

properties of the Magnetic Cataclysmic Variables (MCVs) can generally be explained 
by a model where the field is at most only partially restructured due to accretion. If 
there are systems that have reached the bottom field, they may be found among the 
dwarf novae, and be expected to exhibit quasi periodic oscillations. 

Key words: stars: novae, cataclysmic variables; stars: dwarf novae; stars: binaries: 
close; stars: white dwarfs; stars: magnetic fields. 



1 INTRODUCTION 

In the Magnetic Cataclysmic Variables (MCVs) which com- 
prise some 25% of the Cataclysmic Variables (CVs), there 
is direct or indirect evidence for a magnetic field that in- 
fluences the accretion flow on to the white dwarf. MCVs 
tend to fall into two distinct classes: the AM Herculis-type 
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variables (or polars) and the intermediate polars (IPs) (see 
Warner 1995 for a comprehensive review; Patterson 1994; 
Wood et al. 2000; Wickramasinghe and Ferrario 2000, here- 
after WFOO). 

In the vast majority of the polars, the magnetic field of 
the white dwarf is strong enough to magnetically lock it into 
synchronous rotation with the orbital period through mag- 
netic interaction with the secondary star. Thus, in polars the 
spin period Ps equals the orbital period Porb . Only a handful 
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Figure 1. Magnetic field distribution of the polars (solid his- 
togram), and the intermediate polars (shaded histogram). The 
IP magnetic fields are calculated by assuming that these systems 
are at their equilibrium spin period (Norton et al. 2004), so the 
fields are not measured directly like in polars. The dashed his- 
togram shows the distribution of fields in the high field isolated 
magnetic white dwarfs. 



of polars are slightly asynchronous {Ps/Poih ~ 0.99), likely 
due to recent nova events during which synchronism was 
broken. In polars, accretion occurs directly on to the mag- 
netic white dwarf (MWD) without the formation of an accre- 
tion disc. The magnetic fields in polars are well determined 
through the observation of cyclotron and Zeeman features, 
and are found to be in the range 7 - 230 MG (WFOO). 

The IPs are characterised by the presence of asyn- 
chronously rotating white dwarfs with Ps/Porh ~ 0.01 — 0.9. 
In these systems, accretion occurs via a truncated accre- 
tion disc and magnetically confined accretion curtains (e.g. 
Ferrario & Wickramasinghe 1993; Ferrario, Wickramasinghe 
& King 1993). A direct field measurement is only available 
for one IP; V405 Aurigae for which low-resolution circular 
spectropolarimetry has revealed the presence of cyclotron 
harmonics corresponding to a field of 31.5 MG (Piirola et 
al. 2008) which is in the range of fields found in the po- 
lars. However, as a class, IPs tend to be at the lower end 
(~ 5 — 20 MG) of the field distribution in polars, as it is in- 
ferred by the circular polarization survey of IPs conducted 
by Butters et al. (2009). Fields estimated on the assumption 
that the white dwarf is in spin equilibrium with an accretion 
disc generally confirms this expectation (e.g. Norton, Wynn 
& Somerscales 2004) . The field distributions of the polar and 
the IPs with fields determined either directly (in the case of 
polars), or with the assumption of spin equilibrium and/or 
broadband circular polarization measurements (in the case 
of IPs) are shown in Figure [T] 

There are some parallels and differences with the ob- 
served properties of the isolated Magnetic White Dwarfs 
(MWDs) which fall into two disjoint groups: the high field 
magnetic white dwarfs (HFMWDs) with fields in the range 
10^ — 10^ G, and the low field magnetic white dwarfs 
(LFMWDs) with fields < 10^ G (see WFOO and Wickra- 
masinghe & Ferrario 2009 for a comprehensive review) . The 
HFMWDs comprise some 12 - 15% of all WDs detected in 
volume-limited surveys, while the LFMWDs form a group 
that is at least as numerous. The field distribution of the 
isolated HFMWDs is compared with that of the MCVs and 
is also shown in Figure [T] Although it is this group that is 
most likely to be associated with the MCVs, we note that 
there are distinct differences in the two distributions. The 
peak of the distribution for isolated HFMWDs occurs at 



~ X 10^ G compared to 3 x 10*^ for the MCVs. Furthermore, 
very high magnetic field systems (10* — 10®) are not as well 
represented among the MCVs. 

The magnetic properties of the isolated LFMWDs are 
less well established due to current limitations of measur- 
ing fields ~ 10^ G. It is possible that all WDs are magnetic 
at some level and thus belong to the LFMWD group or it 
could be that most stars are essentially non magnetic. The 
LFMWDs are likely to be associated with the vast major- 
ity of the CVs that shows no direct evidence for coherent 
pulsations at the spin period of the white dwarf as expected 
if there is no magnetically channeled accretion on to the 
surface of the white dwarf. On the other hand, the study 
of Quasi-Periodic Oscillations (QPO) has indicated the pos- 
sible presence of weak magnetic fields in some dwarf no- 
vae (e.g. Warner & Woudt 2002, 2008; Warner & Pretorius 
2008). 

The differences between the observed field distributions 
of the HFMWDs and MCVs may be due to observational 
selection effects or be caused by accretion which affects the 
surface magnetic field of the white dwarf. Isolated HFMWDs 
show no evidence of field decay on the time scale of the 
age of the galactic disc and these conclusions are generally 
consistent with theoretical estimates of time scales for free 
Ohmic decay of the dipolar mode. On the other hand. Gum- 
ming (2002, 2003, 2004) calculated the Ohmic decay modes 
for spherically accreting white dwarfs with a liquid interior 
allowing for compressional heating, and found Ohmic decay 
time scales of several billion years, comparable to the orbital 
evolutionary time scales of typical CVs. These calculations 
appear to suggest that field decay may play a role in our un- 
derstanding of the observed properties of MCVs. However, 
it is unclear whether the assumption of spherical accretion 
is applicable to polars and IPs where accretion occurs via 
funnels onto magnetic polar caps. 

In this paper, we ask the question whether the observa- 
tions of IPs and Polars can be used to distinguish between 
the following three possibilities; (i) accretion induces a a net 
decrease in polar field strength and changes the field 
structure, (ii) accretion screens the field during phases of 
accretion but the field subsequently re-emerges to almost the 
original value when accretion ceases, or (iii) that accretion 
has little or no effect on the field strength and structure. Our 
study is based on a model of non-spherical accretion on to an 
initially dipolar magnetic field which was first developed for 
accreting neutron stars. Here, we extend these calculations 
to conditions appropriate to accreting white dwarfs. The 
model predicts that the magnetic polar cap widens as ma- 
terial is accreted with the field being advected towards the 
equatorial regions by the ensuing hydromagnetic flow (e.g., 
Romani 1990; Konar & Choudhury 2004; Payne & Melatos 
2004; Zhang & Kojima 2006, hereafter ZK06). In the ide- 
alised case where Ohmic diffusion or MHD instabilities are 
assumed unimportant, the polar fleld reaches an asymptotic 
value (the so-called "bottom fleld") as the equatorial fleld 
becomes fully submerged by the advected matter. One of the 
interpretations of the observed properties of neutron stars 
in the Low Mass X-ray Binaries (LMXBs) and the binary 
millisecond pulsars (BMSPs) is that such a bottom field has 
indeed been reached at a value of ~ 10* G (Wijnands & van 
der Klis 1998). 

The paper is organized as follows. We describe the ba- 
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sic model in section 12.11 where we present and discuss the 
equations and solutions for the evolution of the polar field of 
the white dwarf in the case where diffusive effects or MHD 
instabilities are considered unimportant. The expected time 
scales for Ohmic diffusion and advection for conditions ap- 
propriate to accreting white dwarfs are discussed in section 
12.21 The application to the MCVs is presented in section [3] 
where we discuss the various accretion regimes and the pos- 
sible evidence for field restructuring and decline. Finally, 
our main conclusions are summarized in section |4] 



2 THE BASIC FIELD ADVECTION MODEL 

2.1 Advection of polar magnetic zone 

ZK06 developed a simple analytical model which encapsu- 
lates the basic physics of non-spherical accretion onto a mag- 
netised neutron star with a centred dipolar field. On this 
model, as matter accretes along dipolar field lines onto the 
polar cap regions, the accreted matter moves both down- 
wards into the deeper layers of the star, and in a perpendic- 
ular direction dragging the field lines with the flow towards 
the equatorial regions. The net effect is an increase in the 
width of the polar cap region which yields a decrease in the 
polar field strength B. 

Following ZK06, we assume that as accretion proceeds, 
and the polar field B decreases, the magnetospheric radius 
-Rm continues to be given by the standard expression (Ghosh 
& Lamb 1979) for a dipolar magnetic field 



(t>RA = 1.2 X lQ'^{cm)<f)M-^'^ b1^'^ R, 



2/7R4/7pl2/7^-l/7 



(1) 



where Ra is the spherical Alfven radius, and is a pa- 
rameter which is estimated to be ~ 0.5 for disc accretion, 
but is model dependent (e.g. Li & Wang 1999). Further, 
m = M/Mq is the white dwarf mass M in units of solar 
masses, R^ is the white dwarf radius in units of lO" cm, 
Mi6 and Bi are the accretion rate in units of 10^^ g s^^ and 
field in units of 10^ G, respectively. The bottom field B{ is 
reached when the magnetospheric radius Rm ~ R, that is. 



Be = 2.8 X lO'^ (G) Mli^m^'^Rl 



-7/4 



(2) 



In the ZK06 model, the approach to the bottom field 
state is described by assuming that the motion that drags 
and restructures the field occurs in the regions of the star 
that are built up by accreted matter. It is further assumed 
that when the mass of this layer exceeds a critical value 
AMcrit, the accreted matter becomes incorporated into the 
interior of the star, and does not partake in further restruc- 
turing of the field. These assumptions, combined with the 
equations of continuity and magnetic flux conservation, lead 
to the following equation for the polar field B in terms of 
the initial polar field Bo before accretion, and the amount 
of mass AM accreted 



B = 



{l-[C/exp(y)-l]^}V4 
2CAM 



(3) 



, AM = Mt, C = 1 + VI - xo and xo = 



where y 

" 7AMcrit 
/ B \ R 

I -p— J — — — , where Rmo is the initial magnetospheric 
\BoJ Rmo 

radius. The parameter ^ allows for slippage of field lines 



across matter with ^ = 1 corresponding to the ideal case of 
flux freezing. 

Following ZK06, and by analogy with neutron stars, we 
assume that AMcrit is determined by the condition that the 
density at the base of the accreted layer equals the mean den- 
sity of the white dwarf = 2x lO^(M/O.6M0)i?g^ (gs"^). 
Hydrostatic balance gives the pressure beneath a layer of 
mass AM to be 



AM g 

= 2.1 X 10^' (erg cm"^) 



(4) 



— \ R-" 
O.IMq J \Mq " 



where g — GM/R? is the acceleration due to gravity and it 
is assumed that the pressure is that of a degenerate electron 
gas with an equation of state (Shapiro & Teukolsky 1983) 

P = ifp^/^ = 6.8 X 10"^° (erg cm-3)p^/=' , (5) 

where — p/lQ^ {g / cm?). The mass density beneath this 
layer is 



/ AM ff ■y'/^ 
= 2 X 10^ (erg cm" 



AM 



0.1M( 
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M 

M0 
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and its thickness is 
/ AM 



H = 



0.85 X 10* (cm) 



K a"'/" 



(6) 
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(7) 



AM 
0.1 A/q 



2/5 



M 

M© 



-3/5 



R. 
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The average density of a white dwarf of mass M = Q.QMq 
is p = 2 X 10^ g cm~"^ so that with our assumption we have 
AMcrit ~ O.1M0. This attained at a depth H/R ~ 0.1 below 
the white dwarf surface. 



2.2 Advection vs Ohmic diffusion 

In the previous section, we discussed non-spherical accretion 
in the advective limit with the parameter ^ allowing for the 
possibility of inefficient coupling of matter on to field lines in 
a global sense. However, we expect that Ohmic diffusion and 
MHD instabilities will act to counter the advection of the 
field differently in different regions of the star (e.g. equatorial 
and polar) so that the use of a global parameter ^ is too 
simplistic. 

Generally, MHD magnetic field evolution is described 
by the induction equation which includes the contributions 
of two primary terms, MHD flow and magnetic diffusivity. 
The relative importance of these terms can be measured by 
the magnetic Reynolds number 7?.em which is the ratio of the 
Ohmic diffusion time to the flow characteristic time which 
can be estimated from 

tdiff VL 



where V is the flow velocity, L is the length-scale over which 
the field changes, r] = c? /A-kg is the magnetic diffusivity, 
and (J the electrical conductivity. In our case, we can take 
the flow velocity to be the velocity with which the radius 
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of the star shrinks as the mass of the star increases as a 
consequence of the mass-radius relationship of white dwarfs 
(e.g. Shapiro & Teukolsky 1983), 



M 



)Afl6i?9"'P6"' 



(9) 



47r7?2p 

where pa = p / 10^' (g / cm'') . Gumming (2002) has argued 
that accreting WDs should have liquid interiors, with the 
electrical conductivity set by collisions between the degen- 
erate electrons and the non-degenerate ions. For a non- 
relativistic degenerate gas, and for a composition appro- 
priate to CO white dwarfs, approximately (e.g. Yakovlev 
& Urpin 1980; Itoh et al. 1983; Schatz et al. 1999), a = 
6.4 X 10^° (s~^)p6- The magnetic Reynolds number then be- 
comes 



^ aML „ , ,v „_i L 



(10) 



For the polar field, diffusion occurs laterally across the sur- 
face, so an appropriate length scale is L ~ J?. For TZcm 3> 1 
(TZem <C 1), the accretion fiow (Ohmic diffusion) will dom- 
inate the field evolution process. Thus, for accretion rates 
above a critical rate of about 



Mbr(p) = 1.4 X 10 



(11) 



advection will dominate over diffusion and the bottom-field 
state will be reached when a critical mass mass AAfcrit/^ 
is accreted. For significantly lower accretion rates (M <C 
Afbr(p)), the field will diffuse back towards the pole faster 
than it will be advected towards the equator, and accretion 
will have a minimal effect on the polar field strength. 

On the other hand, for the equatorial field it would be 
appropriate to take a length scale L = H, the thickness 
of the region that partakes in the dragging of field lines 
(the electric current zone). With H ~ 0.1-R, we find that 
the equatorial field will be effectively buried above a criti- 
cal accretion rate of A/br(e9) = 1.4 x 10^® (g s~^) or, more 
generally, setting L = H and using equation [T] 



Mbr(eg) = 2.4xl0'''(gs-')i?^/'^ 



M 



3/5 



AM 
0.1 A/q 



-2/5 



(12) 



Above this rate, we expect a bottom field to be reached 
when a critical mass AAfcrit/C is accreted. Somewhat below 
this rate, the field will diffuse outwards faster than it will 
sink into the equatorial region, and we may expect the field 
to be only temporarily enhanced in the equatorial zone. 

It is also instructive to consider separately the time 
scales for accretion and Ohmic decay. The former is given 

by 

(13) 



Af 



6.3 X 10* (yr) 



M 



AM 



O.lAfe J \10iBg s^ 



If we consider diffusion out of a region of length scale L, we 
obtain 



= 3xl0^"(yr)(-j psRt 



For L ^ R, the Ohmic decay time scale of the dipolar com- 
ponent of the field is fdiff (p) ~ 3 x 10^" yr. On the other 
hand, for the equatorial field, L ~ O.li? and the decay time 



scale is ~ 3 X 10* yr. More generally, the Ohmic diffusion 
time-scale across an accreted mass AAf in the equatorial 
zone is 



tdiff (eg) 



4.4 X 10* (yr) 



AM V'^^f M 



Mc- 



0.1 A/, 



-3/5 

R. 



(14) 
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3 APPLICATION TO MAGNETIC 
CATACLYSMIC VARIABLES 

In the strongly advective limit - that is at high enough accre- 
tion rates where Ohmic diffusion cannot counter the advec- 
tion of the field, we may expect the properties of accreting 
white dwarfs to fall into the following broad regimes. For 
these estimates, we assume that initially, the field is strong 
enough to satisfy the condition xq <^ 1 (or R <^ Rmo), and 
also set C = 1. 



of the 



accreted 



• Low values 
(AM < O.OOIMq) 

If AM < AMb < AMcrzt (with xo < 1) ,the solution to 
equation (|3} is given approximately by 



B 



Bo 



(1 + i^)^/- 



(15) 



where AMb = 0.5a;oAArcrit = 0.5(7?/Rmo) AAfcWt. For a 
high field white dwarf in an MCV {Bo ~ 10* G and Bi ~ 
10'^ G), we expect the above approximate solution for AM < 
AMb ~ 10"^AfQ. 

• Intermediate values of the accreted mass 
(AM ~ O.OIM0) 

If the accreted mass satisfies the condition AMb ^ AM <C 
AMcrit, then we have the following approximation from 
equation (|3}- 



B = 0.8B{ 



( AM„ 



7/4 



(16) 



V AM 

Therefore, equation (|16p implies that the infiuence of the 
initial magnetic field on the magnetic evolution has little 
effect at this stage, while the WD magnetic field is scaled 
by its bottom field B{. 

• High values of the accreted mass (AM ~ O.IMq) 
The bottom field Bf is reached when the accreted mass 
approaches AM/ given by 



AM/ = 3.51og(C)AA/crit ^ O.25(Ar0) 



(17) 



\ 0.1 Af0 

Further accretion has no effect on the value of this field. 

To summarise, the results of calculations of the evolu- 
tion of the WD magnetic field with accreted mass are plotted 
in Figure [2] These show that the polar field decreases as the 
accreted mass increases, and saturates once ~ 0.1 A^q has 
been accreted and a bottom field of 3 x 10"^ G has been 
reached. The solution is influenced by the initial field only 
when the total accreted mass is less than ~ O.OIAz/q. Be- 
yond this point, the final (bottom) field depends almost ex- 
clusively on the accreted mass, and the bottom field value 
of 3 X 10^ G is independent of the initial field. 
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Figure 2. The magnetic field versus accreted mass diagram with 
various initial field strengths from Bo = 10^ G to Bo = 10® G as 
labeled in the figure, where the mass accretion rate is set to be 



ZK06 have argued that as the bottom field is reached, 
the stellar field lines are redistributed such that, before be- 
coming buried, the surface equatorial field will approach 
Beq — > B + {R/H){Bo - B) and the polar field will re- 
duce to B ^ Bo yielding Bcq — > {R/H)Bo ~ lOBo if we 
take H ~ O.IR. This means that in an MCV where the field 
has been significantly affected by accretion, a complicated 
field structure may be expected. 



3.1 Is there evidence for field evolution in MCVs? 

In order to assess whether the above solutions could be ap- 
plicable to the MCVs, we must consider the current ideas on 
the orbital evolution of CVs and MCVs. CVs are envisioned 
to evolve to lower orbital periods following first contact (ini- 
tial Roche lobe overflow after common envelope evolution). 
The loss of orbital angular momentum is driven either by 
magnetic braking (MB) by the wind that emanates from the 
secondary, and/or by gravitational radiation (GR). A typical 
system that first comes into contact (say with Porb ~ 8 hr) 
will evolve to Porb ~ 3 hr driven mainly by MB, at which 
point the secondary star becomes fully convective. Accord- 
ing to the canonical model, MB is drastically reduced at 
this point, and the star, which is by now out of thermal 
equilibrium, shrinks within its Roche lobe and mass transfer 
ceases. Q. The orbit continues to shrink due to the loss of an- 
gular momentum by GR. Mass transfer re-commences when 
Porb ~ 2 hr and the system continues to evolve to shorter 
periods until the period minimum is reached (see Warner 
1995 for a review) . Thus, the period range Porb ~ 2 — 3 hr, 
in which there is an apparent lack of CVs, is referred to as 
the "CV period gap" . The evidence for a period gap is not as 
strong in the polars (e.g. Webbink & Wickramasinghe 2002 
and references therein). 

^ The reduction in MB cannot be attributed to the loss of mag- 
netic field as had originally been envisaged, since recent observa- 
tions have shown that fully convective M stars also host large scale 
magnetic fields (Morin et al. 2008). The reasons remain unclear. 



As a first approximation, one can assume that the secu- 
lar evolution of a CV is driven by a mean mass transfer rate 
which is determined by the loss of angular momentum by 
MB and GR. The mass transfer rate will then depend mainly 
on the orbital period of the system. However, there are sig- 
nificant uncertainties in formulating a MB law that is appli- 
cable to secondary stars of the masses found in CVs. The- 
oretical studies of the orbital evolution of CVs have shown 
that mass transfer rates of ~ 7 — 12 x 10^^ (g/s) are required 
at the upper end of the period gap if the observed width of 
the period gap is to be reproduced (e.g. McDermott & Taam 
1989; Hameury, King & Lasota 1991). From the different MB 
laws that have been considered that satisfy this requirement, 
we estimate that a CV that first comes into contact with 
a secondary mass A/2 ~ 0.5Mq and a primary mass 
Ml = O.7M0 (Porb ~ 5 hr), Mmb ~ a few x 10^^ g s^V 
On the other hand, Mgr ~ a few x 10^^ g s^^ for typical 
systems in the period range Porb ~ 3 — 8 hr. Also relevant to 
the following discussion is the time scale to cross the 2 — 3 hr 
period gap without mass transfer which can be estimated to 
be tgap ~ 4 X 10* yr. 

The lack of a pronounced period gap in polars has led to 
the suggestion that these systems evolve much more slowly 
due to a drastic reduction in the rate of MB. The reduc- 
tion is attributed to the trapping of the stellar wind in the 
combined magnetosphere of the two stars in synchronized 
systems (Li & Wickramasinghe 1998). As a consequence, 
the secondary remains in near thermal equlibrium during 
its orbital evolution, and when the field is strong enough to 
severely curtail MB, no period gap is formed (see Webbink 
& Wickramasinghe 2002 for results of detailed calculations). 
As a good first approximation, we can assume that the or- 
bital evolution of most polars is driven by angular momen- 
tum loss at a rate that is closer to the GR rate than to the 
standard MB rate. The mass transfer rate due to GR is typ- 
ically of a few ~ 10^^ g s"'^ at Porb ~ 3 hr which more than 
an order of magnitude smaller than the MB rate at the same 
period. 

A disc system such as an IP, will transfer matter at a 
rate ~ 10^^ g s~^, that is well bf within the regime where 
an advective solution may be appropriate. Such a system 
will accrete a mass AM ~ O.2M0 by the time it reaches the 
period gap, if its initial period of contact is ~ 5 hr and mass 
transfer is conservative. Thus, a bottom-field state could be 
reached. On the other hand, a typical polar will accrete at a 
reduced rate of 10^"" g s~^ and no field evolution or decrease 
in polar field strength will be expected, even though it 
may eventually accumulate a similar amount of mass. 

However, it is by no means clear whether mass will 
accumulate on the white dwarf during the secular evolu- 
tion of a CV at the mass transfer rates estimated above. 
The orbital evolution is expected to be punctuated by nova 
explosions on time scales of ~ lO** — 10^ yr during which 
possibly most of the accreted mass will be ejected from the 
white dwarf (Townsley & Bildsten 2004). If this is the case, 
the bottom-field state may not be achieved even in systems 
that accrete via a disc at a high rate. 

In cases where significant field restructuring has oc- 
curred, we expect to observe MWDs in CVs with enhanced 
equatorial field strengths relative to what is expected for 
a centred dipole configuration. This may be refiected in 
the presence of dominant higher order moments in the ob- 
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Figure 3. Magnetic field versus orbital period diagram for tfie 
polars (filled circles) and IPs (filled triangles). The empty tri- 
angles represent IPs which are very close to synchronism. Note 
that the IP magnetic fields are not directly measured, but simply 
inferred (values from Norton et al. 2004). 



served surface field distributions. Detailed observations 
of white dwarfs in MCVs have enabled field struc- 
tures to be determined only for a handful of systems. 

Although there is evidence for non-dipolar field structures 
in many polars, and for the presence of dominant multipolar 
components in some cases (e.g. Beuermann et al. 2007; Rein- 
sch et al. 2005), such large field variations across the stellar 
surface are clearly not observed. Thus, if accretion-induced 
field restructuring has occurred in polars, it must be at a 
very modest level. Unfortunately, we have no knowledge of 
field structure in disc systems, including the IPs where we 
may expect the largest re-structuring/decline of the polar 
field. 

We show in Figure [3] the observed distribution of mag- 
netic field strengths in MCVs as a function of orbital period. 
If we focus on polars, and assume that they typically evolve 
from long orbital periods through the period gap as MCVs, 
there is no evidence for an evolution from high fields to low 
fields. Indeed, if at all, the evidence appears to be in the 
opposite direction. 

If we look at the entire class of MCVs (IPs and Polars) , 
we note that the ratio of the number of polars to IPs is much 
higher below the gap as compared to above the gap. This has 
often been interpreted as evidence in support of the hypoth- 
esis that IPs above the period gap may evolve into AM Hers 
below the period gap. It is indeed the case that the time- 
scale to cross the period gap is similar to the diffusion time 
scale tdiff (eg) for an accreted mass AM ~ AMcrit ~ 0.2AfQ 
(see section [23J. This suggests that a submerged equatorial 
field in a system that has evolved as an IP above the period 
gap may diffuse outwards during the period gap when accre- 
tion ceases. For AM <C AMcrit, as a first approximation, we 
may simply assume that accretion screens the field (Cum- 
ming 2005). We may then expect the field to re-emerge to 
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Figure 4. The spin period of WD versus accreted mass diagram 
with various initial field conditions, where the mass accretion rate 



near its original value in the period gap. On the other hand, 
when A A/ ~ AA/crit, the screening approximation will no 
longer be appropriate, and we may expect the field to re- 
emerge with a significantly decreased polar field strength. 
The apparent increase in the ratio of AM Hers to IPs as 
one crosses the period gap could partly be due to screening 
above the period gap. However, the dominant effect may be 
the greater ease with which systems can lock into synchro- 
nism below the gap due to the combination of a lower mass 
transfer rates and a shorter orbital period (Wu & Wickra- 
masnghe 1991). 

It is possible that systems that evolve from sufficiently 
long orbital periods may accrete sufficient mass even in the 
presence of nova eruptions to reach a bottom-field state. 
Therefore, this depends on the efficiency with which mass 
is ejected in a nova explosion which remains uncertain. We 
predict that if such systems exist, they will form a sub-class 
of CVs with discs with distinctive homogeneous properties. 
Their polar field strengths would be expected to be Bf ~ 3 x 
10'^ G and they should have strongly non-dipolar magnetic 
field structures. 

Their spin frequency should be a fraction ip of the Kep- 
lerian spin frequency at the magnetospheric radius _Rm ~ R'- 



..^O.Om^){±)A, P^33(s)(^^)(i) 



where 



A = 



M 



1/2 



R 



109 cm 



-3/2 



(18) 



(19) 



Observations of Low Mass X-ray Binaries (LMXBs) have 
shown that the maximum spin frequency of a neutron star 
is 619 Hz, which is roughly half the measured maximum kHz 
QPO frequency 1330 Hz, that is conventionally ascribed to 
the Keplerian frequency (van der Klis 2000; 2006). By anal- 
ogy, we may assume ip ^ 0.5 which suggests that systems 
that reach the bottom-field state should have spin periods 
of ~ 30 s which is close to the fastest spin periods observed 
in CVs. We show in Figure |4] the relation between accreted 
mass and spin period for different values of the initial mag- 
netic field -Bo at a fixed mass accretion rate of lO'^^ g s~^. 



Magnetic fields are expected to play a role in explaining 
QPO behaviour in the accreting neutron stars in Low Mass 
X-ray Binaries (LMXBs) . In many respects, the properties of 
the QPOs observed in the dwarf novae parallel what is seen 
in neutron stars in the LMXBs (see the review by Warner 
& Woudt 2008 and references therein) and so it may be 
tempting to identify the dwarf novae (DN) that exhibit QPO 
oscillations with MCVs that have reached the bottom-field 
state. 

It should be noted that only 25% of CVs are magnetic 
and, if at all, only a small fraction of these systems could 
have reached the bottom-field state. Thus, one should expect 
that the vast majority of the DN are likely to be systems 
that are born with low fields ( ~ 10'' G) and a distribution 
similar to that inferred for the low field isolated magnetic 
white dwarfs (see next section). 

3.2 The field distributions of MCVs and isolated 
magnetic white dwarfs 

As noted in the introduction and illustrated in Figure [1] 
there are distinct differences in the field distribution of 
MCVs and of isolated HFMWDs. The peak of the distri- 
bution for MCVs occurs at 3 x 10^ compared to ~ 10^ G 
for the HFMWDs. Generally, compared to the HFMWDs, 
there is a dearth of low field systems and the very high 
field systems are not as well represented among the MCVs. 
The lack of high field systems may be simply related to 
the extreme reduction that is expected in MB, and hence 
in the mass transfer rate in such systems which will make 
them faint and thus much more difficult to detect (e.g. Li, 
Wu & Wickramasinghe 1994). The apparent lack of low-field 
MCVs, however, is harder to explain. 

In assessing the statistics, we note that there are cur- 
rently about 40% of polars with no field determinations. It 
is possible that observations will reveal that most of these 
systems have magnetic fields between ~ 2 — 20 MG. If we as- 
sume that this is the case, and that the fields are distributed 
as a Gaussian in the logarithm with a mean {log)_B = 6.8 
and a standard deviation aiog b ~ 0.25, we obtain the results 
shown in the bottom panel of Figure [5] The distribution of 
magnetic fields in MCVs and HFMWDs then become sim- 
ilar, and a strong case cannot be made for field decay. On 
this scenario, it has to be assumed that most systems that 
are observed as IPs above the period gap will synchronize in 
the period gap and be seen as polars below the gap thus ex- 
plaining the observed reduction in the ratio of IPs to polars 
below the gap. 

On the other hand, if observations show that the 40% 
of polars with no field determinations have fields similar to 
those measured in the remaining 60% of polars, then the 
reduction in the ratio of IPs to polars below the gap can be 
attributed to effects of accretion on the field. 

The magnetic properties of the isolated LFMWDs and 
their incidence are less well established, due to current lim- 
itations of measuring fields ~ 10^ G. Thus, it is possible that 
either all WDs are magnetic at some level or that most WDs 
are essentially non magnetic. In any case, the LFMWDs are 
likely to be associated with the vast majority of the CVs 
that show no direct evidence for coherent pulsations at the 
spin period of the white dwarf as would be expected if there 
were magnetically channelled accretion on to the surface of 
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Figure 5. Top panel: the magnetic field distribution of the 
HFMWDS. Middle panel: the magnetic field distribution of the 
MCVs. Here the field strength of most IPs has been inferred on 
the assumption that the white dwarf is in spin equilibrium with 
an accretion disc. Bottom panel: the magnetic field distribution 
of the MCVs obtained by assuming that 40% of polars with no 
field determinations have magnetic fields between ~ 2 — 20 MG 
(see text). 

the white dwarf . Unfortunately, we have no direct evidence 
of the field strengths of the white dwarfs in such CVs with 
discs (e.g. in dwarf novae). 



4 SUMMARY AND CONCLUSIONS 

The anomalously low fields that have been deduced for the 
neutron stars in the majority of LMXBs (10* — 10^" G) com- 
pared to the magnetic fields in the isolated radio pulsars 
(10^^ — 10^* G) is seen as strong evidence for accretion- 
induced field decay (e.g. Manchester 2006). In these disc ac- 
cretors, the competition between the advection of the field 
from polar to equatorial regions, and the tendency of the 
field to re-emerge by Ohmic diffusion results in a net reduc- 
tion of the field at least during phases of accretion. On the 
other hand, if the LMXBs do evolve into the BMSPs as is 
generally assumed (but see Hurley et al. 2009; Ferrario & 
Wickramasinghe 2007), one can also conclude that the field 
does not re-emerge when accretion ceases, at least not on 
the time scale of the typical age (~ a fewlO® yr) of a BMSP. 
The latter conclusion would imply that accretion does not 
simply lead to a screening of the magnetic field which sub- 
sequently re-emerges to near its original value, but would 
yield a change in the current forming regions, and thus to 
field decay. 

We have investigated whether a similar process may be 
in operation in accreting MWDs in the CVs. In particular, 
we have asked whether there is evidence for the existence of 
group of systems that have reached a bottom field among 
the CVs. 
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We have argued that the mass transfer rates that de- 
termine the secular evolution of CVs is different for systems 
that evolve as polars or as IPs. Based on our simplified model 
for the evolution of the field due to advection, and estimates 
of times scales of Ohmic diffusion, we conclude as follows. 



ACKNOWLEDGEMENTS 

This research has been supported by NSFC (No. 10773017) 
and National Basic Research Program of China 
(2009CB824800). 



(i) In the advection dominated regime, we expect that 
as matter accretes on to the polar caps, the field lines will 
be dragged from the polar caps towards the equator. The 
bottom-field state is reached when ~ 0.2Mq is accreted and 
the equatorial field becomes buried below the white dwarf 
surface. Systems that enter a bottom-field state will have no 
memory of initial conditions (initial field and spin period) 
and are expected to form a homogeneous group with Bf ~ 
3 X 10^ G, and P ~ 30 s. 

(ii) The advection dominated solution is expected to be 
applicable only for accretion rates above a critical value 
Merit ~ lO^^g/s. Below this rate, the field will diffuse out- 
wards faster than it is advected. The polars, with typical 
secular mass transfer rates of a few times 10^^ g s~^, are 
therefore not expected to exhibit significant field restruc- 
turing and reduction in polar field strength. This is 
consistent with the general lack of evidence for field evolu- 
tion with orbital period in polars. 

(iii) The advection dominated solution is expected to be 
applicable to the IPs above the period gap, which typically 
have accretion rates of ~ 10^ g s~^. Systems that first come 
into contact at PorbS hr will accrete AM0.2Mq as they reach 
the upper end of the period gap (Porb ~ 3 hr) if mass transfer 
is conservative. However, when nova eruptions are taken into 
consideration, the mass that is built up by accretion may be 
significantly smaller. We may therefore expect only a small 
proportion of IPs, if any, to reach the bottom-field state. 

(iv) If there are MCVs that have reached the bottom- 
field state, they are likely to be found among the dwarf 
novae that exhibit QPOs. They may be expected to have 
very peculiar field structures dominated by strong equato- 
rial field enhancements. However, this prediction cannot be 
easily verified, since we have little information on the field 
structures of white dwarfs in disc accreting systems. 

We conclude by noting that there are two major differ- 
ences between LMXBs and CVs that make field restruc- 
turing and /or decay more likely in the LMXBs. Due 
to vastly different magnetic moments of the neutron stars 
(typically lower by a factor ~ 10®), there are no magneti- 
cally phase locked systems similar to the polars. Thus, while 
the majority of the CVs that are recognised to be magnetic 
are the strong field polars, which we have argued are unlikely 
to exhibit accretion induced field restructuring due to 
their intrinsically lower mass accretion rates, all the LMXBs 
(which include both high and low field neutron stars) are 
disc systems with intrinsically higher mass transfer rates, 
and are therefore more likely to exhibit field restructur- 
ing and /or decay. A second major difference is that 
thermonuclear runaways similar to nova explosions 
in CVs are not expected to result in mass loss from 
neutron stars, because of their deep gravitational 
potential, so that the critical mass is likely to be 
achieved more easily during LMXB evolution. 
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